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Geometric phase via adiabatic manipulations of the environment 
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We show that geometric phases may be generated in a quantum system subject to noise by adiabatic 
manipulations of the fluctuating fields, e.g., by variation of the system-environment coupling. For a two-state 
quantum system we express this phase in terms of the geometry of the path, traversed by the slowly varying 
direction and amplitude of the fluctuations. We discuss the origin of this phase and possibilities to separate 
it from the known environment-induced modification of the Berry phase. 
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Introduction. After the discovery of the geometric 
phases in coherent quantum systems [H [2] (see also [3] 
on geometric phases in earlier work), it was natural to 
ask, whether these phases can be observed in quantum 
systems coupled to an environment. In particular, since 
the environment typically has a continuous spectrum, 
the gap in the spectrum of the system+environment 
vanishes, which blocks manipulations at frequencies be- 
low the gap. Hence, this criterion of adiabaticity needed 
to be revisited. This is of special interest in connec- 
tion with recent experimental observations of the Berry 
phase and the effect of the environment in supercon- 
ducting circuits [H [5] . 

Furthermore, the concept of a geometric phase needs 
to be defined (extended) for an open system, where the 
wave function and its evolution phases are ill-defined. 
While various formal generalizations have been dis- 
cussed, here we take an operational approach and ana- 
lyze physically relevant, directly observable quantities. 
Specifically, for an adiabatically manipulated isolated 
quantum system the Berry phase is a phase of a matrix 
element of the evolution operator for the wave function. 
For a system, weakly coupled to noise, we analyze the 
corresponding evolution operator of its density matrix, 
which determines observable quantities. We find that 
this operator contains phase factors, which reduce to 
the standard Berry phase for isolated systems. 

Analysis of this kind allows one to define geometric 
phases for open systems, and they can be compared to 
their values in the limit of the vanishing coupling to the 
environment. An important question is whether and 
how are they modified by the noise. 

We analyze this problem, using Bloch-type master 
equation for the evolution of the density matrix. Let us 
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comment on some earlier attempts to analyze this evo- 
lution operator (cf. the discussion in Ref. 6j). In var- 
ious contexts, Bloch-type master equations have been 
used. However, if the standard master equation is used 
in the case of an arbitrary (slow) time variations of the 
Hamiltonian, one finds the same value of the geometric 
phases as without the noise. We stress, however, that 
the dissipative terms in the standard Bloch equations 
are derived for a static Hamiltonian, and they should 
be modified (rederived), if the Hamiltonian varies in 
time, even slowly. This analysis shows, that the noise- 
induced terms (Tj- and T2-terms as well as the "Lamb 
shift" for a spin system) are modified by the variations 
of the Hamiltonian which leads to modification of the 
geometric phases. 

In Ref. [7] this study has been performed for a spin- 
1/2 system in a magnetic field (any two-level system 
reduces to a spin-1/2), which varies along a cone, with 
the direction of the noisy contribution to the field along 
the axis of the cone. It was found that the Berry phase 
acquires a noise-induced contribution. In Ref. [5] the 
analysis was performed for arbitrary loops, traversed by 
the tip of the magnetic field, and it was found that the 
modification of the Berry phase is of geometric origin 
(similar to the Berry phase itself) and has a quadrupo- 
lar symmetry. Moreover, it was observed that this mod- 
ification of the Berry phase is complex, which implies a 
geometric contribution to dephasing. 

Here we describe another possibility to generate an 
environment-induced geometric phase. The analysis in 
Refs. [5] was done under the assumption that the 
spin-bath coupling is fixed. However, typically adia- 
batic manipulations of the Hamiltonian are performed 
by changing control parameters of the quantum system, 
and this can easily influence the coupling to the environ- 
ment or the properties of the bath (cf. the effect of flux 
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noise in Ref. [5]). The strength and the matrix struc- 
ture of the system-bath coupling may also be modified 
deliberately. Under these circumstances new contribu- 
tions to the Berry phase may arise, and the analysis of 
earlier work cannot be applied directly. Here we ana- 
lyze the geometric phases in situations, when not only 
the Hamiltonian but also (or only) the coupling to the 
environment is varied adiabatically, and find the corre- 
sponding contributions to the geometric phase. 

Hamiltonian. A two-level system is equivalent to a 
spin-1/2, and we use the spin notations to describe its 
dynamics. The Hamiltonian of a spin coupled to a fluc- 
tuating field can be written as 

W = -~B<T--i"n<T + W elaV) (f) 

where the field B corresponds to the controlled part of 
the Hamiltonian. The second term describes the in- 
fluence of the environment, the fluctuating field. To 
demonstrate the effect and following earlier work we 
consider uni-directional fluctuations, but allow the di- 
rection to vary slowly in time. Thus X is the bath 
operator and the slow-varying n(t) indicates the direc- 
tion and the strength of the fluctuations. The last term 
in Eq. (fTJ) describes the dynamics of the environment 
and, in particular, determines the correlation functions 
of X. We assume that the field X has zero average (in 
the absence of coupling to the spin) . 

Variations of n may result from the dynamics of the 
environment or changes in the spin-environment cou- 
pling, ft may be induced deliberately or be a side effect 
of the adiabatic variations of the spin Hamiltonian B. 
To illustrate the effect, here we analyze the slow varia- 
tions of the direction and the strength of a unidirectional 
noise with fixed correlations; this analysis can be easily 
generalized to include more general fluctuating fields, in 
particular, with a varying power spectrum 2 ). 

We analyze the phase accumulated by the system 
(and the dephasing) between times and tp. Detec- 
tion of this phase may involve preparation of the initial 
state (e.g., a superposition of f) and J.}) and the final 
direct or indirect measurement. We do not specify de- 
tails of these events (see below). We further assume 
the following conditions for the time scales involved: 
t c ,_B _1 <C ip,T2, where r c is the noise correlation time 
and T-i is the dephasing time scale. This implies, in par- 
ticular, that the noise is weak and short-correlated [9], 
and that on the time scale tp of the evolution the noise 
correlations are local. Note also that the coherence de- 
cays on the time scale T2, and at longer times tp ^> T2 

2 'The results dlOI I. mil l below apply also in this case, with the 
substitution S(fi) -> S{Q,t) = f dre iUT (i [X(t + f), X(t- §)]+). 



the phase information is exponentially suppressed. B 
and n vary at typical frequencies ~ lu — 2n/tp, and the 
adiabaticity parameter is lo/B. 

Reference frame. To find the effect of the time vari- 
ations of the Hamiltonian B and the noise n, it is con- 
venient to make a time-dependent transformation of the 
wave functions or, in the spin language, a transforma- 
tion to a rotating reference frame (RF), in which the di- 
rection of the B-field is stationary p] [6] . In that frame 
we use the standard procedure to write the master equa- 
tion for the density matrix and thus find the evolution 
operator. 

Below we find that the evolution of the off-diagonal 
entry of the density matrix in the RF, p^, is described 
by the expression 

PTlW-PTllOK^^W^ ^- (2) 

Here the real part of the phase gives the angle of the spin 
precession about B, while its imaginary part (due to T 
and 5$) describes the decay of the transverse spin (de- 
phasing) . In Eq. ([2|) the first term in the exponent gives 
the dynamical phase and dephasing; this term scales 
with the total time tp. The second term is the geometric 
phase, $0 + 8<f>, insensitive to time reparametcrization: 
the bare Berry phase $0 depends only on the geometry 
of the path B(i), and the environment-induced complex 
contribution depends on the geometry of B(i) and 
n(t). " 

The analogy with the spin dynamics shows that, sim- 
ilar to the Berry phase for a closed system, the (total 
and) geometric phase at time tp for an open system is 
meaningful only when the direction of the field B at 
tp coincides with its initial value (however, there is no 
restriction on n(tp)). At the same time, the imaginary 
part of the phase, the decay of coherence, is well-defined 
also for open paths (cf. Ref. [6]). 

To choose the frame, we choose eigenstates |ft), ||t) 
of B(t)cr for each t (we omit the subindex t below). 
In the spin language this fixes the axes of the RF: 
z = B/£>, while x is such that |f t } + || t ) is the +1- 
eigenstate of xcr. We have the freedom to choose the 
phases of the states \\ / J.), or the direction of the x, y- 
axes; this choice does not affect the results, we only 
assume that the states (the axes) vary slowly enough, 
at frequencies <C B; we further suppose that the x,y- 
axes assume their initial values if we consider a closed 
loop, i.e., when the direction of B(i) at the final time 
returns to its initial value. 

In this rotating frame the pseudo-magnetic field is 
B' = B + uj, where u> is the angular velocity of the RF. 
The direction of B' is stationary to the leading order 
in the adiabatic parameter but differs from the z-axis, 
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and it is convenient to choose it as the z'-axis in the 
RF ( u x' j/'z'-frame" ) . This transformation of the basis 
in the RF, which changes the direction of the third axis 
from B to B', is only weakly time-dependent (uj ~ u 2 ), 
and thus corrections to the pseudo-magnetic field at this 
step are negligible. 

The eigenstates of the Hamiltonian in the RF (that 
is, the eigenstates of B'<r = B'cr z /) can be easily found 
as 



IT') = lt) + *^ll) 
II') =U)-i^lt) 



(3) 



(notice the choice of the overall phases). 

Using these eigenstates we find the level splitting in 
the RF: 

S'«S + i((TrT>-U|i>) • (4) 



To proceed with the calculation in the RF, we need to 
find the components of the fluctuating field. We intro- 
duce the coordinates n z , n± 
in the xyz-fr&me: n<x = irt+(7_ 



n x ± in y of the n- vector 
^n_(7+ +n z a z , where 



a± = a x ± ia y . One finds that its coordinates in the 
primed frame are slightly different, and to the leading 
order in u>: 



n+ 



2in z 
~B~ 



(I It), 



n„ - n z = 



B 



n+(T ||) +n_U If) , 



(5) 
(6) 



and n'_ = (n' + )*. 

Below we find that under the specified conditions the 
equation of motion for the off-diagonal entry p-^ of the 
density matrix in the primed basis is decoupled from the 
other matrix elements, and the evolution factor contains 
information about the Berry phase [10] . This factor can 
be extracted from the measurements of the spin state. 

Evolution of the density matrix. To calculate Pn(t), 
we derive a Bloch-Redfield master equation for the evo- 
lution of the reduced density matrix of the two-level 
system considered (cf., e.g., Refs. [TT] [12j [T3J [9] and ref- 
erences therein). Let us outline the derivation: we be- 
gin from the Liouville equation for the density matrix 
of the combined system and environment, assuming an 
initially factorized density matrix, then make the trans- 
formation to the interaction representation, and expand 
the evolution operator to the second order in the system- 
bath coupling. The derivation is performed under the 
assumption of weak and short-correlated noise, which 
allows one to reduce the integro-differential equation of 
motion to a markovian differential equation. Finally, 
the secular (or rotating-wave) approximation allows us 



to decouple the evolution of different entries of the den- 
sity matrix. 

After the expansion in the perturbation V' = 
— ^Xn'a-' and averaging over the state of the environ- 
ment, wc find: 



(id t + B') Pu (t) 



, , <[[|| t )(T t |,vmv'(*i)])*i. 

-DO 



(7) 



(We can set the lower limit of integration to — oo for 
t ^> t c , noise correlation time, which is also the con- 
vergence scale of the integral in Eq. ([7]) [Til H21 H3] ; the 
~ r c -intervals at the boundary contribute to the bound- 
ary phase <5$b, see below.) 

Using the secular approximation (r <C B), we ob- 
tain: 



(id t + B')p u = -ip n / S(t-h) x 

J —OO 

(^(ty+ftiJe'^i^^ + niWn^tx)) dh , (8) 



where S(t-ti) = ±(X(t)X(ti)+X(ti)X(t)) is the sym- 
metrized correlation function of the noise. 

For an isolated system, i.e., neglecting the rhs of 
Eq. JSJ), one finds the phase acquired as $total(£p) = 
B'dt = f* P (B + (uB)/B)dt = $l P Bdt + <P° BP , 
which is the sum of the dynamical phase and the con- 
ventional Berry phase. -$g P (mod 2ir) is given by the 
solid angle subtended by the path B(f). 

The fluctuations on the rhs of Eq. © give rise to 
the dephasing of the element and the noise-induced 
contribution to the phase. If we neglect all order-B, h 
effects, this expression produces terms oc t in the total 
phase: setting n'_(t) = n-(t) = ri-{t\), n',(t) = n+(t), 
and B' = B, we find: 



r = i 



dn s(n) 

2^n- B + iO 



n!/^L,. (9 ) 



2vr tt + iO 



Here S(£l) is the Fourier transform of S(t). Note that 
r is complex, its real part gives the decoherence rate, 
while its imaginary part determines the modification of 
the level splitting ("Lamb shift"), cf. Ref. [9]. 

To find the noise-induced modification of the Berry 
phase, we expand the rhs of Eq. ([5]) to the first order in 
B and h. Integrating the evolution over time, we find a 
geometric contribution: 
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<5$ = 



dft S(Q)(3B - 2n) \ nB n(B x dB 



4tt B(Q - B + iO) 2 J B 
dil S(fl) \ B(n x dn) 
B 



B 2 



2ir (Q - B + iff) 2 
dG + (5$ b ■ 
Here the quantity G is given by 



(10) 



G(n,fl) = G(|n H 

>+| 2 /"dn 



\,n z ,B) = 

s(fi) 



(ii) 



2tt (n- B + iO) 2 



dn s(Q) 
2^ (n + ioy 



The last term <5<I>b is the boundary contribution and will 
be discussed below. Eq. (|10[) is our main result. 

Let us discuss various terms in Eq. ([TO]). The first 
term gives the contribution due to variation of the 
Hamiltonian B at fixed n, which was found in Ref. [BJ. 
Here we present it in a form, independent of a choice of 
coordinates, which may be convenient for applications. 
We see from this expression that this contribution arises, 
when the vector of the pseudo-magnetic field B rotates 
around n. 

Geometric phase generated by the rotation of noise. 
The next term of Eq. IjlOp. J . . . dn, describes the cor- 
rection to the geometric phase due to the variations in 
the system-environment coupling. It appears when the 
direction n of the fluctuating field rotates around the 
bare field B. We refer to this term as the term due to 
rotation of the noise. 

Let us illustrate the origin of this term by consid- 
ering the simplest example, in which this contribution 
arises: the bare Hamiltonian retains its initial value, 
B = const, whereas the direction n of noise varies in the 
plane, transverse to B, with angular velocity p (FigH^). 
Then according to Eq. (fT0[) the first term vanishes. We 
postpone the discussion of the last line in this equation 
and consider the contribution due to rotation of noise. 
The considerations are simplified by going to the frame, 
rotating together with n around B. In that frame n 
is constant; the field B is also constant but it acquires 
an extra contribution, B — > B + p. This implies that 
the rate T @ is modified compared to its value in the 
laboratory frame, and hence the total phase acquires a 
contribution i J STdt = i J(dB^)pdt: 



,*.noisc 
gcom 



l dB dLp - 



(12) 



The phase <l?g°om nas the following geometric inter- 
pretation: for a closed loop n(i) and B = const this 
phase is given by the flux of the uniform field (see 

Fignw 



b n = e 2 



dn 



s{n) 



2tt (CI- B + iO) 2 



(13) 



through the surface spanned by the loop n(i). This re- 
sult follows directly from the Eq. (|10[) . 

These results follow immediately from the well- 
known expression of the Berry phase via the solid an- 
gle, if the noise is slow (typical 1] C B). Indeed, 
during the rotation, transverse adiabatic fluctuations 
of the total magnetic field B + Xn sweep the area 
i|n_L| 2 (X 2 )p = ±\n + \ 2 (X 2 )p per unit time (see FigQJi). 
Thus the geometric phase is: 



4> 



slow noise 



l 

= — \n>- 
2 



! (X 2 



dp 
B 2 ' 



(14) 



in agreement with Eq. (|10p . 

Thus, we find that the variation of the system-bath 
coupling may generate a geometric phase, that is, a 
phase which depends only on the geometry of the path 
n(t). This phase has both real and imaginary parts, that 
is, it describes a geometric contribution to dephasing as 
well. 

One can also understand the contributions in 
Eq. (fTO)) in the following way: if neither B nor n rotate 
about each other, i.e., when the (B, n)-plane retains its 
direction, the contribution to the geometric phase (due 
to variations of B and n within the plane, that is, due to 
changes of their magnitudes and the relative angle) can 
only be of potential nature [6J and give the term J dG 
in Eq. (fT0|) . The generic case can always be reduced to 
it by going to a proper rotating frame; for instance, one 
can begin with the RF used in the calculations above 
and add an additional rotation about B' to keep the 
phase of n + constant. This additional rotation, with 
the angular velocity tp (here n+ = |rt+|e l(p ), modifies 
the field as B' — » B' + tp and thus the dephasing 
rate in Eq. ([9]), producing a contribution to the phase 
5$ <— i J dpds^, which coincides with the second term 
of Eq. (jTOj (cf. also QUI). 

The boundary phase. The last term in Eq. (| 10|) is a 
boundary contribution and is accumulated in the vicin- 
ity of the initial and the final points, t = and t = tp; it 
cannot be presented as a line integral. Let us comment 
on the origin and the value of this term. In the deriva- 
tion of the master equation ([7|) we assumed a factorized 
density matrix of the system + environment at an ear- 
lier time to = ~ oo, the latter served as the lower limit 
of integration in Eq. ([7]). In fact, the exact initial condi- 
tions for the density matrix and, in particular, the exact 
value of to do not matter except for times t very close 
to to, in an interval ~ r c . However, the behavior of the 
density matrix and the phase accumulated within this 
short time interval are sensitive to the initial conditions; 
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this contribution to the phase depends on the details of 
the preparation at the beginning of the adiabatic ma- 
nipulations and does not scale with the total time tp 
of the Berry-phase experiment, hence this boundary ef- 
fect should be formally ascribed to the geometric phase. 
The same considerations apply to the short final inter- 
val, when the (quantum) measurement of the final state 
is performed, the exact value of the phase being sen- 
sitive to the details of the read-out procedure. These 
two contributions form the boundary term ^$b- For 
an abrupt measurement at t = tp the final boundary 
term vanishes; the initial boundary contribution also 
vanishes for the slow or preliminary preparation, which 
corresponds to to = — oo; in contrast, for an abrupt 
preparation with to — one finds from Eq. ((7]) that 
£$ b = 2G(n(0),B(0)). 

Notice that the boundary phase masks the other con- 
tributions in Eq. ([7|). For a closed path, as emphasized 
in Ref. [6] , the contribution of the first terms can be en- 
hanced relative to the boundary term by traversing the 
path several times, N: their contribution scales cx N, 
whereas <5$b retains its value. However, the term J dG 
vanishes for a closed loop and cannot be enhanced by 
this method. In other words, to calculate and compare 
to experiment the last line in Eq. 1(7} one needs to take 
into account the details of the preparation and read-out. 

In conclusion, we found that slow variation of the 
properties of the noise or of the system-environment 
coupling may result in a geometric phase. We found this 
contribution for a two-level system in a fluctuating field 
with a slowly varying direction. We thank A. Shnirman 
for discussions. This work was partially supported by 
the projects INTAS 05-1000008-7923, MD-4092.2007.2, 
and the Dynasty foundation. 




FIG[T] a) The rotation of noise about a constant mag- 
netic field B generates a geometric phase given by the 
flux of a uniform field b n through the loop n(t). b) The 
area swept by the vector B + n(t)X(t) in the plane, per- 
pendicular to B 



9. Yu. Makhlin, G. Schon, and A. Shnirman, in: New Di- 
rections in Mesoscopic Physics (Towards Nanoscience), 
eds. R. Fazio, V. F. Gantmakher, and Y. Imry, pp. 197- 
224, Springer, 2003. 

10. Our results (S.V. Syzranov and Yu. Makhlin, in: Elec- 
tron transport in nanosystems, eds. J. Bonca and S. Kru- 
chinin, Springer, to appear) show that in multi-level sys- 
tems the situation may be more complicated. 

11. F. Bloch, Phys. Rev. 105, 1206 (1957). 

12. A.G. Redfield, IBM J. Res. Dev. 1, 19 (1957). 

13. C.W. Gardiner and P. Zoller, Quantum Noise, Springer, 
2000. 



1. M.V. Berry, Proc. R. Soc. Lond. 392, 45 (1984). 

2. A. Shapere and F. Wilczek, eds. Geometric phases in 
physics, World Scientific, 1989. 

3. G.B. Malykin and S.A. Kharlamov, Usp. Fiz. Nauk 173, 
985 (2003) [Phys. Usp. 46, 957 (2003)]. 

4. P.J. Leek, J.M. Fink, A. Blais, R. Bianchetti, 
M. Goppl, J.M. Gambetta, D.I. Schuster, L. Frunzio, 
R.J. Schoelkopf, and A. Wallraff, Science 318, 1889 
(2007). 

5. M. Mottonen, J.J. Vartiainen, and J. P. Pekola, e-print 
larXiv:0710.5623l (2007) . 

6. R.S. Whitney, Yu. Makhlin, A. Shnirman, and Y. Gefen, 
Phys. Rev. Lett. 94, 070407 (2005). 

7. R.S. Whitney and Y. Gefen, Phys. Rev. Lett. 90, 190402 
(2003). 

8. G. Falci, R. Fazio, G.M. Palma, J. Siewert, and V. Ve- 
dral, Nature 407, 355 (2000). 



